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We study the kinetics of helix-coil transition in water solutions of gelatin (collagen protein) by 
optical rotation techniques combined with thermal characterization. By examining the rates of 
secondary helix folding, and covering a very wide range of solution concentrations, we are able to 
identify a universal exponential dependence of folding rate on concentration and quench temperature. 
We demonstrate a new concentration-temperature superposition of data at all temperatures and 
concentrations, and build the corresponding master curve. The results support the concept of a 
diffuse helix-coil transition. We find no concentration dependance of the normalized rate constant, 
suggesting first order (single) kinetics of secondary helix folding dominate in the early stages of 
renaturation. 
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Kinetics of protein folding is one of the cornerstone 
problems in understanding the full mystery of biologi- 
cally active macromolecules Q, |^ . A number of exper- 
imental techniques focus on a great variety of natural 
and synthetic polypeptides undergoing their globular col- 
lapse, or denaturation into the coil. Although one might 
want to think of such a transition in terms of individual 
macromolecules, there is much evidence that cooperative 
effects between different chains play a significant role, 
an effect often referred to as oligomerization 3]. This 
extreme variability has, to date, defied any attempt on 
universal description, beyond the classical Zimm-Bragg 
abstraction A\. 

In this Letter we examine perhaps the most well- 
studied protein, collagen, and discover that the complex- 
ity of its folding kinetics, strongly dependent on temper- 
ature and concentration in the solvent, can be dramat- 
ically simplified. All of this kinetics can be scaled onto 
a single master curve by a new procedure we call the 
"concentration-temperature superposition" . There is a 
remarkable analogy: the classical time-temperature su- 
perposition 1^ has allowed master curves to describe 
the glass transition in a variety of thermal viscoelastic 
systems and served a great purpose in rheology for the 
last 50 years. A much more recent discovery of time- 
concentration superposition 0,0 has allowed the univer- 
sal description of dynamic glass, or jamming transition, 
in lyotropic systems such as colloid suspensions. In both 
cases the cooperativity of interparticle interaction is the 
key, and we shall argue that it is relevant for the folding 
kinetics as well. 

We must emphasize that we study the principal helix- 
coil transition in collagen chains, i.e. the folding of sec- 
ondary helices. This is achieved by our experimental 
methods focusing on optical rotation of linearly polarized 
light at a wavelength away from any absorbtion band. It 
has been recently demonstrated that the dominant con- 
tribution to the measured rotation rate arises from the 
secondary helices |9|. This is different from the more 
common studies of gelation kinetics in gelatin, which 
is controlled by the tertiary triple helices and is usu- 



ally studied by rheological methods. The field of sol-gel 
transition and mechanical response of networks is broad 
and well-established [lol [lll | . Recent important observa- 
tions of universal master curves, describing specifi- 
cally the storage modulus in gelatin, also make connec- 
tions with glassy dynamics. We shall see from the results 
below that our present work is not related to these ideas: 
we examine much shorter times at which the secondary 
structure is formed, while the kinetics of subsequent rhe- 
ological is a much slower process. 

Collagen is the main protein component of white fi- 
brous connective tissues such as skin, tendons and bones. 
The fundamental unit of the native collagen structure is a 
rod consisting of three individual molecular strands, each 
twisted into a secondary left-handed helix, which analo- 
gous and very similar to the classical a-helix in other 
polypeptides. In collagen, three segments of secondary 
helix wrap together to form a tertiary right-handed su- 
perhelix stabilized by further interchain hydrogen bond- 
ing . Collagen is extracted from tissues by hydrolytic 
degradation, which denatures it producing the resulting 
material commonly known as gelatin. 

Gelatin is dissolved in water by heating the solution to 
40° C. Above this temperature collagen chains have ran- 
dom coil configuration. On cooling transparent gels form 
containing extended physical crosslinks. X-ray diffrac- 
tion and transmission electron microscopy measurements 
suggest that the crosslinks are formed by partial rever- 
sion to ordered triple-helical segments, separated along 
the chain contour by pep tide residues still in the random 
coil configuration |14Lll5l |. At very low concentrations the 
renaturation, or folding process is completely intramolec- 
ular and proceeds by a back folding of the single chains 
[l6j |. With increasing concentration the renaturation be- 
comes increasingly intermolecular. Gelation occurs at 
concentrations above 5 mg/ml due to the formation of an 
infinite elastic network in the gelatin solution, see ^3,01 
and references therein. 

Gelatin owes many of its uses to this coil-helix tran- 
sition. The kinetics of the transition has been exten- 
sively studied for many years using a variety of tech- 
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niques such as differential scanning calorimetry (DSC), 
scattering and rheometry, as well as optical rotation. Fig- 
ure Gla) shows a sequence of DSC traces of helix-coil 
transition that we obtained in gelatin solutions of differ- 
ent concentration, while Fig.^b) shows the correspond- 
ing transition line, Tm{c). A very flat dependence of 
the transition temperature on concentration in the ther- 
mal phase diagram, extrapolating to c — > 0, indicates 
that the transition is largely intrachain. Note that the 
gelation phase diagram, obtained by rheological methods, 
would show a drop of the gel point Tg{c) below 5 mg/ml. 
Also note that, since the melting temperature depends 
strongly on the gelation temperature and time, a broad 
range of Tm(c) values are found in the literature. 

Gelatin, as any protein, is an optically active material 
in both the random coil and helical states. However, due 
to coherent long-range chiral ordering, helical domains 
rotate the plane of light polarization much more strongly 
than the individual chiral aminoacids in the coil state. 
Thus, the coherent optical activity gives a direct indica- 
tion of the fraction of the monomers in helical states. 

In many biopolymers, including polysaccharides and 
DNA, the coil-helix transition is very fast and resembles 
a true first order phase transition '20'|. In gelatin, how- 
ever, the helix nucleation step lasts several minutes and 
subsequent growth of the helices proceeds even slower, 
at a logarithmic rate [23|. This allows one to study its 
detailed kinetics. Most investigations have focused on 
isothermal renaturation rather than temperature scan- 
ning studies, because of the difficulty in accessing the 
equilibrium state. This slowness of the process will allow 
us to scan and superpose the renaturation kinetics over 
a range of temperatures and concentrations. 

In one of the earlier studies, Flory and Weaver stud- 
ied helix growth rates in very dilute gelatin solutions 
(c < 4 mg/ml) and found first order kinetics J2J, i.e. that 
the folding rate is concentration- independent. They pos- 
tulated the coil-helix transition proceeds via an interme- 
diate state formed by intramolecular rearrangement of a 
single chain. Assuming this state consists of a secondary 
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FIG. 1: (a) DSC traces demonstrating the helix-coil tran- 
sition in gelatin/water solutions of diflterent concentrations, 
labelled on plot, (b) A thermal phase diagram of this tran- 
sition, Tm{c), obtained from the DSC data (0), with results 
for different preparation methods from the literature (□ - fl^ . 
■ - figl'). The value Tc is obtained from the concentration- 
temperature superposition analysis below. 



helix segment, consideration of the minimum stable seg- 
ment length leads to the Flory- Weaver expression for the 
rate constant for the renaturation after quenching the 
dilute solution: 



ki = B exp 



\kTAT J ' 



(1) 



where A and B are constants, T is the quench temper- 
ature, and AT — Tj^ — T is the degree of supercooling 
below the equilibrium melting temperature T^. 

The first order kinetic analysis agrees with other ex- 
perimental observations of dilute gelatin solutions jl6l |. 
In semidilute solutions concentration-dependant kinetics 
has been observed 0, . Consequently, new mecha- 
nisms have been proposed, involving either different seg- 
ments of the same chain or up to three different peptide 
chains having to interact to form stable helices. The 
topological consequences of these arrangements on helix 
formation and gelation are very significant. 

Goddard et al. 18], suggesting an analogy between re- 
naturation and crystallization, modelled the kinetics in 
terms of Avrami exponents |2J|: x = 1 — exp[— fcii"], 
where x is the helix fraction. They obtained the expo- 
nent n close to unity, suggesting one-dimensional growth 
from predetermined nuclei and supporting the first order 
kinetics ideas. Early optical rotation studies confirmed 
that simple exponential kinetic occurs in dilute solutions, 
but only in the initial stages of renaturation )i25t .2dj . 
Most recently, Guo et al. 17] studied initial renaturation 
rates of gelatin solutions at semidilute concentrations up 
to 0.12 g/ml and observed what appeared as a combina- 
tion of first order and second order kinetics, i.e. the rate 
of growth of the normalized helix fraction had a linear 
concentration dependence. This seems unusual, in view 
of the triple nature of tertiary helix linkages in collagen. 
However, the observations of second-order gelation kinet- 
ics, in particular the elastic modulus proportional to c^, 
have been made for over a century ^Tf\ . A two-step mech- 
anism with the rate limiting step formation of a nucleus 
of two helices wrapped together, followed by rapid subse- 
quent wrapping of another coil segment to form the triple 
helix has bee n propo sed a long time ago to account for 
these findings [l^ilSg]- However, we shall argue that there 
is a big difference between optical rotation and rheolog- 
ical studies. The latter return the viscoelastic response 
of the gel and thus rely on the triple-helix linkages of 
the network. Our opinion |^ is that the properly mea- 
sured optical rotation signal is determined by secondary 
helices and not the tertiary structure. We attribute the 
second-order kinetics of Guo et al. 0| to the range of 
examined concentrations - only up to 0.1 g/ml (which 
is also the range in the rheological work of Normand et 
al. ^3)- Overall, examining the literature of the last 30 
years, one finds no universality, nor agreement between 
different concepts. 

In an attempt to link together all of the past findings, 
we studied the helix-coil transition in gelatin over a much 
wider concentration range, up to 0.4 g/ml. The universal 
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FIG. 2: Typical traces of normalized optical rotation, di- 
rectly proportional to the secondary helix density, as func- 
tion of time after quenching the denatured collagen solution 
to Tq = 11°C. Different data sets correspond to concen- 
trations labelled on graph. Solid/dashed lines are fits with 
^(l_exp[-i/T]). 



master curves we construct from the new superposition 
procedure give explicit predictions for c w 0.65 g/ml, and 
can be further extended to the full range of temperatures 
and concentrations. The high-sensitivity differential op- 
tical rotation detector we used to measure the secondary 
helix content has been described in, e.g., [2^. The raw 
measurement of the total angle of polarization plane ro- 
tation, is divided by the (constant) sample thickness 
to produce the rotation rate. It is then normalized by 
the solution concentration to obtain a specific rotation 
[a] — {l/c)[d'^/dz], from which we subtract a bare value 
[a]o corresponding to the average aminoacid optical ac- 
tivity, separately measured in the coil state at high tem- 
peratures. The resulting difference is proportional to the 
concentration of correlated helices in the medium. The 
typical reading, for a range of concentrations, is shown 
in Fig. 121 for a solution quenched from 50°C to the value 
T = 11°C in this case. As a result of [a] normalization, 
the y-axis in Fig. [3 is directly proportional to the helix 
fraction in the sample. 

There are important and delicate issues of the slow 
drift of collagen towards its natural state [23|, reflected 
in the small deviation of the long-time data from the sim- 
ple exponential fit of each data set (shown by the lines 
in Fig. O. For instance, the rheological study of Nor- 
mand et al. [l^ is heavily based on this regime, where 
the tertiary structure (and the gel elasticity) are being 
consolidated and the elastic modulus increases dramati- 
cally. We, however, are concerned with the initial rates 
of secondary helix growth, i?o(c, T), essentially the slopes 
i?o = d[a]/dt = A/t a.t t ^ Q, from the fitted func- 
tions A[l — exp(— t/r)] in Fig. [3 For each quench tem- 
perature, these slopes depend on the solution concen- 
tration. The summary of this dependence is given in 
Fig. 13 The highly non-linear concentration dependence 
of growth rates is apparent from the exponential fits of 
all data sets. We note that the highest solution concen- 



FIG. 3: Initial rates of helical nuclei growth, _Ro at t ~+ 0, 
as functions of solution concentration, for a range of quench 
temperatures labeled on plot. Solid/dashed lines are fits with 
exponential B exp[yc]. The arrow indicates the concentra- 
tion range of earlier studies [l^. 
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FIG. 4: The master curve of growth rate vs. concentration, 
Ro{c), superposing the results for different quench tempera- 
tures (labelled on plot). The arrows point at different data 
sets; the solid line is a universal fit with B exp[Y'c|). 



tration quantitatively studied so far was c ~ 0.12 g/ml, 
by Guo et al.\vi\. 

Guided by the experience in time-temperature and 
time-concentration superposition, we notice that the sets 
of data for the growth rates in Fig.|21can be shifted along 
the concentration axis by an amount that is a function 
of quench temperature. Selecting a reference tempera- 
ture (at this stage, we arbitrarily choose Tj-d — 21°C) 
one scales the concentration c for each data set such that 
c = (3 ■ c, with the coefficient (the shift factor) a function 
of the quench temperature for each data set, /3 = (3{T). 
The fact that these sets do superpose means that there 
is a universal underlying expression for the growth rate 
i?o[/3(T)c]. The resulting master curve, Fig.01 is an indi- 
cation of such a universal physical process that controls 
the helical nucleation and growth in all regions of the 
(T, c) -phase diagram. 

It is important to establish a law of the shift factors 
[3 dependence on temperature. Figure [3 gives the plot 
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FIG. 5: Dependence of the shift factor /3 on quench temper- 
ature Tq. The reference temperature chosen at 21°C corre- 
sponds to the original, non-shifted data - or /3 = 1. The linear 
fit (dashed line) crosses the /3 = line at ~ 30° C. 



of a set P{T) required to produce the superposed mas- 
ter curve for the initial growth rate. It follows quite an 
obvious linear function, /? = 0.1 (Tc — T), with the "crit- 
ical temperature" Tc ~ 30°C (more precisely, 302. 5°K) 
defined as the point at which (3 = 0. This temperature 
is labelled in the gelatin phase diagram and the sequence 
of DSC scans, Fig.QJa). Clearly, its relation to the tran- 
sition line is not accidental. In fact we should conclude 
that (within errors of our experiment and analysis) this 
critical temperature is indeed the line of the helix-coil 
transition, Tc = Tni(c). With this assumption, the he- 
lical growth rate, as obtained from the master curve in 
Fig. O follows the equation 



d[a] 



= B cxp [5(T„, - T) c] 



(2) 



with fixed parameters B « 0.06deg/dm (g/ml) ^s^^ 
and h « 0.6 (ml/g)/°K. Note that we find no critical or 
singular behavior near the transition. The growth rate 
becomes exponentially small at T > Tm. In retrospect, it 
is perhaps not surprising that we did not obtain a critical 
vanishing of the growth rate, as would be the case with 
classical phase transitions, or was suggested by the old 
Eq. Q . Our present finding supports the idea of a diffuse 
nature of helix-coil transformation in a single chain, with 
chain fluctuations capable of creating a non-vanishing 
helical fraction above Ti„(c) and disrupt thermodynami- 
cally equilibrium helices formed below this temperature. 
In other words, there is no sharp phase boundary Ti„(c) 
for the helix-coil transition. In equilibrium models of this 
transition one identifies T,„ as a point of steepest gradient 
in helical fraction x{T), while in our study of kinetics we 
obtain roughly the same value at a point when Rq{c,T) 
changes from the decaying to the growing exponential 
function exp[l'c]. 

Having established the universal master curve for the 
folding rates i?o(c, T), we must now make contact with 
the other kinetic studies, that traditionally focus on the 
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FIG. 6: Normalized folding rates k as functions of concentra- 
tion at different quench temperatures (labelled on plot) . 



normalized helical fraction, in our notation defined as 



xit) 



No 



A(c,T) 



1 — exp(— i/r) 



(3) 



according to Fig. |21 and earlier work The origin 

of the exponential concentration dependence, Eq. jSJ, 
is in the constant A{c), while the normalized rate k = 
1/t = Rq/A. The concentration dependence of this rate 
is a signature of process topology, e.g, a constant value 
means first-order (single-molecule) folding process. Guo 
et al. claimed a second-order process, i.e. a linear 
dependence k = ki + k2C + but we clearly arrive at 
different conclusions. Figure shows our results for the 
normalized rates k over a concentration range much wider 
than in any previous work. There are some delicate issues 
of fitting the data that includes the slow long-time tail, 
which certainly accounts for the large noise in the k{c) 
data. However, the qualitative picture clearly suggests 
that there is no relevant concentration dependence and, 
therefore, the single-chain folding process (first-order) ki 
is in fact dominant at early stages of renaturation. This 
may appear a surprising result as, with increasing concen- 
tration, intramolecular rearrangements of single chains 
become limited. Also, the triple-helical nature of colla- 
gen helices is well established indeed 0|. However, let 
us remember that our experimental method inherently 
looks at the secondary helices (the signal from triple he- 
lical structures is very low due to the large wavelength 
mismatch). The observed first-order (i.e. single-chain) 
nature of the principal folding process appears to be an 
indication that the transition from coil to secondary he- 
lix is a single chain process, while the association of these 
helices to form the triple network linkages occurs subse- 
quently and contributes less to the helical onset than was 
previously thought. 

In conclusion, by applying a new concentration- 
temperature superposition we found a universal master 
curve describing the initial folding rates Rq of gelatin so- 
lutions over a broad range of concentrations and temper- 
atures, spanning the nominal helix-coil transition point 
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Tm- The exponential (c, T)-dcpcndcncc of i?o arises from 
the increasing overall amount of helices in different sys- 
tems. In contrast, the normalized helical fraction 
grows with a c-independent rate, suggesting the first- 
order kinetics. One needs to study other oligomerizing 



proteins to test the further universality of the discovered 
master curves. 
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